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ABSTRACT 

We propose a possible way to detect baryons at low redshifts from the analy- 
sis of X-ray absorption spectra of bright AGN pairs. A simple semi- analytical 
model to simulate the spectra is presented. We model the diffuse warm-hot in- 
tergalactic medium (WHIM) component, responsible for the X-ray absorption, 
using inputs from high-resolution hydro-dynamical simulations and analytical 
prescriptions. We show that the number of OVII absorbers per unit redshift 
with column density larger than 10^^'^ cm^^ - corresponding to an equivalent 
width of ~ 1 km/s - which will be possibly detectable by XEUS, is ^ 30 per 
unit redshift. Constellation- X will detect ~ 6 OVII absorptions per unit red- 
shift with an equivalent width of 10 km/s. Our results show that, in a ACDM 
Universe, the characteristic size of these absorbers atz~0.1is^l/i~^ Mpc. 
The filamentary structure of WHIM can be probed by finding coincident ab- 
sorption lines in the spectra of background AGN pairs. We estimate that at 
least 20 AGN pairs at separation 20 arcmin are needed to detect this fila- 
mentary structure at a Scr level. Assuming observations of distant sources using 
XEUS for exposure times of 500 ksec, we find that the minimum source flux 
to probe the filamentary structure is 2 x 10~^^ erg cm~^ s^^, in the 0.1-2.4 
keV energy band. Thus, most pairs of these extragalactic X-ray bright sources 
have already been identified in the ROSAT X\\S\iy Survey. Re-observation of 
these objects by future missions could be a powerful way to search for baryons 
in the low redshift Universe. 

Key words: Cosmology: theory - intergalactic medium - large-scale structure 
of universe - quasars: absorption lines 

1 INTRODUCTION 

The census of baryons in the low redshift Universe shows that a large fraction of them has not been detected yet. 
While almost 80% of the baryons reside in the Lya forest at 2 ~ 2 (Ranch 1998), estimates at lower redshift imply 
that the observed baryon density is significantly lower (Fukugita et al. 1998; Nicastro et al. 2002b) and does not match 
the nucleosynthesis constraints (Buries & Tytler 1998). Hydro-dynamical simulations suggest that a large fraction of 
baryons at z ~ could be in a warm-hot intergalactic medium (WHIM) , with a temperature between 10^ and 10^ K 
and with moderate overdensities 5 ~ 10 — 100 (Ostriker & Cen 1996; Cen & Ostriker 1999; Dave et al. 2001; Cen et 
al. 2001). Detection of this component is particularly difficult and requires very sensitive UV and X-ray satellites. 
Recently, various attempts to model analytically the WHIM component have been made. Perna & Loeb (1998) 
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used the Press & Schechter theory (Press & Schechter 1974) and an isothermal model for the gas component to 
address the detectability in absorption of the hot and highly ionized intergalactic medium (IGM), which resides in 
the outskirts of clusters and groups of galaxies. Valagoas, Schacffcr & Silk (2002) showed that the tcmporaturc-dcnsity 
relation found in hydro-dynamical simulations can be reproduced using straightforward physical arguments. In this 
simple picture the IGM consists of two phases: the cool photo-ionized component which gives rise to the Lya forest 
and the WHIM, dominated by shock heating. 

Hydro-dynamical simulations of structure formation in the context of cold dark matter models have been the 
major contribution to the theoretical understanding of the WHIM. Hellsten et al. (1998), referred to the network 
of filaments and sheet-like structures, possibly soon in the spectra of background bright sources in the X-ray band, 
as the 'X-ray forest', in analogy with the Lya forest. Their main result is that OVII and OVIII absorptions by 
the intervening IGM can be seen by future X-ray missions. Cen et al. (2001) used large box-size hydro-dynamical 
simulations to detect OVI absorption by the WHIM in the UV band and find results in reasonably good agreement 
with observations. Fang et al. (2002a) performed hydro-dynamical simulations to study the metal distribution in the 
IGM and, assuming collisionally ionized gas, compared their results with semi-analytical models based on the Press 
& Schechter formalism. Kravtsov et al. 2002 addressed WHIM detectability using hydro-dynamical simulations to 
model the gas distribution in the Local Supercluster region and imposing constraints from the MARK HI catalogue 
of galaxy peculiar velocity. 

More recently, Chen et al. (2002) have made an extensive study of the detectability of OVII and OVIII absorbers 
using the z = output of a hydro-dynamical simulation of a ACDM Universe, and assuming a uniform metallicity of 
0.1 solar for the IGM. Their results confirm that the detectability of these absorbers is challenging. It has been shown 
that the number of intervening OVI absorbers per unit redshift is quite high at low redshift (Tripp et al. 2000). Thus, 
a very promising technique to detect OVII and OVIII systems is to search for corresponding UV OVI absorptions in 
the spectra: this allows to lower the S/N ratio needed for OVII and/or OVIII detectability. Nicastro et al. (2002a) 
presented the X-ray 5a level detection by Chandra of a resonant absorption from WHIM along the lino of sight 
towards the blazar PKS 2155-304. Moreover, Fang et al. (2002b) reported a detection of an OVlll absorption line 
along the sight-line towards PKS 2155-304 with Chandra. They constrained the gas that gives rise to this line to 
have an overdensity in the range 30-350 and a temperature of 4-5 x 10^ K, in agreement with hydro-simulations. Also 
emission by the WHIM could be important and several models have been proposed (e.g. Scharf et al. 2000; Phillips 
et al. 2001; Zappacosta et al. 2002). Here, we will focus on absorption: the bulk of the missing baryons should reside 
in a medium which is too rarefied to be detectable through its X-ray emission or Sunyaev-Zel'dovich effect (Ostriker 
& Cen 1996). 

These very recent observations are very encouraging and show that we have just started to detect the first 'trees' 
of the X-ray forest. Future X-ray missions such as XEUS and Constellation- X will be much more accurate a probe of 
the current cosmological models and will allow a detailed study of the WHIM properties. 

In this paper we will present a semi-analytical technique to simulate X-ray absorption spectra. There are two 
advantages in using a semi-analytical technique instead of the more accurate hydro-dynamical simulations. First, 
semi-analytical models are not computationally expensive and are not limited by resolution and box-size effects. 
Second, semi-analytical techniques can be very efficient when coupled with detailed hydro-dynamical simulations 
when the latter can be used to calibrate the former and allow one to explore very quickly the parameter space. 
For example, it will be straightforward to study the effects of the metallicity of the IGM, different UV and X-ray 
backgrounds, scatter in the temperature-density and in the metallicity-density relation, peculiar velocities etc. For 
the estimates wo want to make here, semi-analytical techniques, which use the inputs of hydro-dynamical simulations, 
will show to be successful in reproducing potentially observable quantities such as the column densities of OVII and 
OVIII absorbers. 

In particular, we will focus on the simulations of absorption spectra of AGN pairs and we will try to determine 
the characteristic size of the absorbers by identifying absorption features in the simulated spectra. The estimate of 
the extent of the absorbers, using the information contained in the transverse direction, could be a very powerful 
probe of the filamentary network of WHIM. An extensive analysis of the correlations of Lya clouds aX z > 2 using 
hydro-dynamical simulations has been made in Miralda-Escude et al. (1996) and Cen & Simcoe (1997). McGill (1990), 
Crotts et al. (1994), Bechtold et al. (1994), Charlton et al. (1997) performed similar analysis for observed Lya clouds 
and showed that the coherence length in the transverse direction of the clouds could be significantly larger than 
the size of a single cloud. This argument has been clarified recently using simple physical arguments (Schaye 2001) 
and hydro-dynamical simulations (e.g. Theuns et al. 1998), which indeed show that the clouds are embedded in a 
network of filaments. The idea is very simple: while a single linc-of-sight (LOS) allows to probe only the size of a 
single absorber, pairs of LOS are sensitive to the random orientation of the filaments and diffuse blobs and can better 
sample the WHIM. Another advantage of using the information contained in a sample of pairs is that this method 
minimizes the contamination of the WHIM signal we have to detect with the intrinsic lines of the background source. 
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This is due to the fact that the positions of the intrinsic absorption or emission features of the source along the LOS 
are random. Thus, their effect can be statistically modelled and removed as it has been done recently to recover the 
linear dark- matter power spectrum from simulated Lya absorptions (Vicl et al. 2002a). 

The outline of the paper is as follows. In Section 2 we describe the ingredients needed to simulate the X-ray 
AGN spectra focussing both on the inputs from hydro-dynamical simulations (subsection 2.1) and on the analytical 
modelling (subsection 2.2). In Section 3 we test how well the semi-analytical model matches the properties of simulated 
lines extracted directly from hydro-simulations. This comparison is done in terms of the predicted number counts of 
absorbers per unit redshift range. In Section 4 we discuss a possible way to recover the characteristic size of these 
absorbers by using the so called 'hits-and-missos' statistics and discuss detectability of these absorption features by 
satellites of the next generation such as XEUS and Constellation- X. Section 5 contains our main conclusions. 



2 SEMI-ANALYTICAL MODEL OF THE X-RAY FOREST 

This Section presents the semi-analytical method used to simulate X-ray spectra of AGN. In Subsection 2.1 we will 
present the hydro-dynamical simulation used, focussing on the basic ingredients which will constitute the inputs for 
the semi-analytical model: i) the probability distribution function (pdf) of the gas; ii) the gas mctallicity - density 
relation; Hi) the gas temperature - density relation; iv) the X-ray background. In Subsection 2.2 we will describe in 
detail the technique used, with the relative mathematical implementation (see also Appendix Al), to produce the 
simulated spectra of AGN pairs. 

2.1 Model parameters from hydro-dynamical simulations 

The hydro-dynamical simulation used here is a ACDM model with the following parameters: Ha = lOO/i km s~^Mpc~^ 
with h = 0.67, Q.om = 0.30, fiob = 0.035, A = 0.70, as = 0.90, and the spectral index of the primordial power spectrum 
n = 1. The box-size of the simulation is 25 comoving Mpc on a uniform mesh with 768^ cells. The comoving 
cell size is 32.6 kpc and the mass of each dark matter particle is ~ 2 x lO'^ Mq (further details can bo found 
in Cen et al. 2001). The simulation includes galaxy and star formation, energy feedback from supernova explosions, 
ionization radiation from massive stars and metal recycling due to SNe/galactic winds. Metals are ejected into the 
local gas colls where stellar particles are located using a yield Y — 0.02 and are followed as a separate variable 
adopting the standard solar composition (some slices through the simulation box at various redshift outputs can be 
found at http://astro.Princeton.EDU/~cen/PROJECTS/p2/p2.html). 

To properly simulate X-ray absorptions we have to know the actual X-ray and UV background of the simulation at 
z = 0, which is computed self-consistently given the sinks and the sources in the simulation box. We have checked that 
there are very small differences between the hydro-simulation background and the estimate from ShuU et al. (1999): 
luv = Iuv{E/lS.6eVy^-^, with I'^y = 2.3 x 10"^^ erg cm"^ Hz"^ sr"^ s"^ for the UV background (which includes 
a contribution from AGNs and starburst galaxies), and Ix = Ix{E/Ex)~^'^''^ exp {-E/Ex), with = 1.75 x 10"^'^ 
erg cm"^ Hz"^ sr"^ s"^ and Ex = 40 keV (Boldt 1987; Fabian & Barcons 1992; Hellsten et al. 1998; Chen et al. 
2002), for the X-ray background. The difference between the amplitude of the UV and X-ray background of the 
hydro-simulations and the estimates above is less than a factor of 2, in the range 0.003 ^ E (keV) ^ 4. In Figure 1 
we plot the X-ray and UV background of the hydro-simulations at 2: = (dashed line) , the estimates of ShuU et al. 
(1999) (long dashed line) and Fabian & Barcons 1992 (dot-dashed line), and the sum of these two (continuous line). 

Among the two backgrounds, the X-ray one plays a major role in determining the statistics of the simulated 
spectra. We will assume that the amplitude of the ionizing backgrounds scales like (1 -|- zf , in the redshift range 
< 2: < 1. This assumption will not influence significantly our results: given the fact that we will generate mock- 
spectra at a < 0.3, this will determine a scaling of the amplitude by a factor of less than 2. 

Along many randomly chosen sight-lines parallel to one of the axes of the simulation box, we compute the gas 
density and peculiar velocity, the gas temperature, and the metal density. This allows us to obtain a sample of 1000 
simulated absorption spectra. 

From the 2 = output we extract the pdf of the gas density, which is shown in Figure 2, this will be implemented 
in the semi-analytical (hereafter SA) model. Starting from the linear density field, generated with the right correlation 
properties predicted by linear theory (Subsection 2.2), we assign each pixel a new non-linear density in such a way 
that the resulting non-linear pdf is that extracted from hydro-dynamical simulations. 

The gas density-temperature relation is shown in panel (a) of Figure 3, where we report the mean value in each 
density bin and the scatter. This result is in reasonable agreement with the recent results obtained by Yoshida et 
al. (2002 - their Figure 2) and also by Dave et al. (1999), who studied the gas cooling processes in the framework 
of Smoothed Particle Hydro-dynamics (SPH) simulations. We note that the scatter in temperature can be very high 
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Figure 1. The X-ray and UV background. The short-dashed line shows the background extracted from the 2 = output of 
the hydro-dynamical simulations and used to compute our spectra. The estimates of ShuU et al. (1999) for the UV background 
are represented by the long-dashed line. The results of Fabian & Barcons (1992) for the X-ray background are also reported 
(dot-dashed line). The sum of the latter two is represented by the continuous line. 
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Figure 2. Probability distribution function of the gas overdensity for the 2 = output of the ACDM model. 

for values of log(l + S) > —1, spanning more than one order of magnitude. We report the gas density-metallicity 
relation in panel (b). The trend here is as expected: low density regions are the most metal poor, while at large 
densities the metallicity can be solar. In this plot the scatter is extremely large for the intermediate density regime, 
i.e. —1 < log(l + S) < 0.5, while it is significantly smaller for larger overdensities. This plot has to be compared with 
similar plots shown in Cen & Ostriker (1999b), but we find that here the scatter is somewhat larger than in previous 
simulations. 

These curves, with the relative scatter, assumed to be drawn from a Gaussian distribution, will be implemented 
into the SA model. 
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Figure 3. Panel (a) shows the temperature density relation for the gas component, while panel (&) shows the metallicity as 
a function of density. The continuous lines represent the mean value in each bin, the dashed lines represent the scatter {la 
values). The quantities shown are volume weighted. 



2.2 Simulation of absorption features in X-ray spectra of AGN pairs 

We start here from the model introduced by Bi and collaborators (Bi et al. 1992, 1995; Bi 1993; Bi & Davidsen 1997 

- hereafter, BD97) for generating a Lya absorption spectrum along a line-of-sight. This simple model predicts many 
properties of the absorption lines, including the column density distribution and the distribution of line widths (6 
parameters), which can be directly compared with observations. 

The BD97 model is based on the assumption that the low-column density Lya forest is produced by smooth 
fluctuations in the intergalactic medium which arise as a result of gravitational growth of perturbations. Linear 
density perturbations in the intergalactic medium Sl'^^{'x.,z) are related to the underlying linear dark-matter (DM) 
perturbations by a convolution, which models the effects of gas pressure. In Fourier space one has: 5}!^^ (k, z) = 
WiGuik, z)D+(z)dQ^'^(k), where -D+(z) is the growing mode of density perturbations (normalised so that -D+(0) = 1) 
and (5j"^(k) is the Fourier transformed DM linear over-density at z = 0. A commonly adopted low-pass filter is 

WiGM{k, z) = (1 + k'^/kj)~^, which depends on the comoving Jeans wave-length \j{z) = 
where fcs is Boltzmann's constant, Tm the gas temperature, /i the molecular weight and 7 the ratio of specific heats, 
fiom is the present-day matter density parameter. BD97 adopt a simple lognormal transformation to obtain the IGM 
density in the mildly non-linear regime from the linear one. 

Several approximations are involved in this modelling. A first point which can bo important is that the Jeans 
length depends on density and temperature, and therefore Jeans smoothing should be an adaptive smoothing of the 
density field, while the Fourier implementation cannot take into account this local aspect (Viel et al. 2002b). In 
addition, a filtering of the dark-matter density field at the Jeans scale may be not very accurate: different filtering 
functions result in very different gas density fields and the real 'filtering scale' which reproduces better the gas 
distribution, in the linear regime, depends on the thermal history of the IGM (Hui & Gnedin 1998; Matarrese & 
Mohayaco 2002; Gnedin et al. 2002). 

The modelling of the gas component using the lognormal transformation and a smoothing at the Jeans length 
would determine a poor description of the gas distribution at 2: ~ 0, where the density fluctuations are in the non- 
linear regime. Even the tight relation between density and temperature for the gas responsible of Lya absorption is 
no longer valid for the warm-hot intergalactic medium. In this case, hydro-simulations suggest that the scatter in the 
T - 5 relation, for the WHIM, can be very large (Dave et al. 1999; Yoshida et al. 2002). 

At least three inputs, taken from hydro-dynamical simulations, can be used to improve the model: i) the 
probability distribution function (pdf) of the gas density; it) the gas density-temperature relation; Hi) the gas 
density-metallicity relation. Viel et al. (2002b) showed that by using directly the pdf of the gas density extracted 
from the hydro-dynamical simulations, the pdf of the flux predicted by the semi-analytical model is in significant 
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better agreement with that in the hydro-simulations. The other two inputs will allow to take into account the physical 
state (temperature and metallicity) of the structures responsible for the absorptions. 

Here wc will follow closely the proscriptions used in Vicl et al. (2002a, 2002b). The modelling is based on the 
generation in Fourier space of ID density and velocity fields, with the right correlation properties predicted by the 
theory. 

Starting from a linear power spectrum P{k) we generate ID LOS linear random fields for the IGM density and 
peculiar velocity, with the low-pass filter window function Wigm- If we draw a LOS in the x\\ direction at a fixed 
coordinate x±_ the Fourier component at redshift z of the linear density contrast and of the linear peculiar velocity 
field of the intergalactic medium (IGM) along the LOS are: 

5jf^(fe||,^|xx) = D+{z) j ^e*— WGM(^fc| + fei,^)5o°''(fe||,kx), (1) 
t;f^(fc||,^|xx) = ihE+{z) j l^e*— ^^^WGM(y'fc^ + fci,^)5o°^(fc||,kx), (2) 

with E+(z) — H(z)f(Q,m,^K)D+{z)/{l + z). Here f{Q.m;^K) = — din D+(2;)/dln(l + z)), and Vim and f^A are the 
matter and vacuum-energy contribution to the cosmic density (e.g. Lahav et al. 1991), respectively; H{z) is the 
Hubble parameter at redshift z: H{z) = i?o-\/f^Om(l + z)^ + f2or-(l + 2)^ + ^oa, with Qor = 1 — ^om — ^oa- In our 
case we use the ACDM parameters of the previous subsection. 

We then simulate the corresponding density and peculiar velocity at a distance r± from the first LOS. The cor- 
relation properties between these 4 fields are computed using an algebraic implementation, as described in Appendix 
Al. 

The linear density field is mapped into the non-linear one using directly the pdf of the gas distribution obtained 
from hydro-simulations. This has been done with a rank-ordering technique between the linear density field and the 
non-linear one, which intrinsically consists in a monotoiiic and deterministic mapping of each generated field. To 
account for the scatter, we assign each pixel a new temperature and metallicity value by adding to the mean value 
(continuous lines of Figure 3) a random number taken from a Gaussian distribution with the same standard deviation 
as the hydro-simulations. 

Peculiar velocities are approximated by those predicted by linear theory, as the peculiar velocity field is known to 
keep linear even on scales where the density contrast gets mildly non-linear. This simplification should not be critical 
in simulating the spectra, as the effect of peculiar velocities results in a shift of the line and in an alteration of the 
profile (BD97), which will have small impact on the recovered column density of the absorber. 

After having simulated the gas distribution, the photoionization code CLOUDY (Ferland et al. 1998) is used to 
compute the ionization states of metals. In this analysis, we have considered both the collisional ionization and the 
photo-ionization. The ionizing background consists of the sum of a UV background, arising from AGNs and galax:ies, 
and a X-ray background, probably produced by a population of moderately X-ray luminous AGN at ~ 2 (Fabian & 
Barcons 1992), while for the soft X-ray part (< 1 keV) a significant contribution comes from the WHIM itself (Cen 
et al. 1995; Phillips et al. 2001). 

Among the heavy elements that can produce absorption lines in the X-rays Oxygen is the most abundant one 
and produces the strongest lines. A more complete treatment should also include other elements such as C, N, Ne, Fe 
and Si but, in our simplified scheme, we will concentrate only on Oxygen lines. For a full list of resonant absorption 
lines we refer to Verner et al. (1996). We simulate absorption spectra for the two strongest transitions of the ions 
OVII {E = 0.57 keV, A = 21.6A, / = 0.7) and OVIII {E = 0.65 keV, A = 18.97A, / = 0.42), with / the oscillator 
strength. 

The density of each ion is obtained through: ni{x) — nnix^XiYz^Z j Zq^ where Xi is the ionization fraction of 
the ion as determined by CLOUDY and depends on gas temperature, gas density and ionizing background, uh is 
the density of hydrogen atoms, Z is the metallicity of the element and Fzg is the solar abundance of the element. In 
Figure 4 we plot the ionization fraction for OVII and OVIII as a function of the temperature for two different values of 
the IGM density corresponding to an overdensity &igm ~ 6 (left panel) and Sigm ~ 120 (right panel), corresponding 
to densities of 1.4 x 10~® cm~^ and 3 x 10~^ cm~^. To check if the metals considered here are in photoionization 
equilibrium we compute the recombination time tree ~ (iiigm Qrec)^, with Orec the recombination rate. For OVII 
and OVIII, in the temperature ranges considered here, one has that a^ec ~ 2 x lO"^"' cm^ s~^ (Mazzotta et al. 1998), 
where both the dielectronic and radiative recombinations are taken into account. This means that the Hubble time 
is larger than the recombination time roughly for densities ^ 10~^ cm~^. Thus, for the WHIM, the approximation 
of photoionization equilibrium should be reasonable. 

Given the density of a given ion along the LOS, the optical depth in redshift-space at velocity u (in km s~^) is 
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Figure 4. Logarithm of the ionization fraction as a function of gas temperature for OVII (continuous line) and OVIII (dashed 
line). The left panel shows the results obtained with CLOUDY for an IGM overdensity 5 = 6, while in the right panel the 

results are for an overdensity S = 120. 



T/H = y dyni{y)V[u-y-vl{y),b{y)] (3) 

where ao,i is the cross-section for the resonant absorption and depends on Aj and fi, y is the real-space coordinate 
(in km s~^), V is the standard Voigt profile normalised in real-space, b = (2kBT/mic^y^^ is the thermal width 
and we assume that vi = vigm- Velocity v and redshift z are related through dA/A = dv/c, where A = Aj(l + 2). 
For the low column-density systems considered here, the Voigt profile is well approximated by a Gaussian: V = 
{y/7vb)~^ exp[— (-u — y — W|j (y))^/&^]- The X-ray optical depth t will be the sum of the source continuum, which we 
assume we can determine, and that of eq. (3). Finally, the transmitted flux is simply = exp(— t). 

We generate 400 spectra with our semi-analytical technique at a median redshift of z ~ 0.05 and ~ 12000 km/s 
long. Both for the semi-analytical model and for the hydro-simulations spectra have been produced with a resolution 
of ~ 3 km s"'^. The goal is to reproduce number and characteristic size of these absorbers with the simulation of 
AGN spectra by a simple semi-analytical recipe, and address their detectability. Figure 5 shows the OVII distribution 
at 2 = extracted from the hydro-simulation. We can see that the OVII connects high density regions and displays 
a network of filaments and diffuse blobs. Prom this Figure we can argue that strong absorptions arise in virialized 
halos, while the weaker ones, possibly dotoctod by future X-ray missions, sample filaments connecting high density 
regions. The OVIII distribution is quite similar to the OVII one (Chen et al. 2002). 



3 TESTING SEMI-ANALYTICAL MODELS WITH HYDRO-DYNAMICAL SIMULATIONS 

In this Section we show how well the semi-analytical model works, after having calibrated it with the inputs obtained 
from hydro-dynamical simulations: the pdf of the gas, the gas density-temperature relation and the gas density- 
metallicity relation (with the relative scatter) and the X-ray background. We compare the column density distribution 
function for OVII and OVIII predicted by the semi-analytical model with those extracted directly from the hydro- 
dynamical simulations. To this aim, the SA method has been run cutting off the large scale power for scales larger than 
the box-size of the hydro simulation {25h~^ Mpc). This allows a more direct comparison between the hydro-simulation 
and the semi-analytical model. We checked that the statistics considered here are only marginally influenced by this 
assumption, by running another set of SA spectra without the cut in the power-spectrum. We conclude that the 
box-size of the hydro-simulations is large enough to study reliably these absorbers. Column densities have been 
computed by integrating directly the simulated optical depth and not using any fitting routine available, both for 
the hydro-simulations and for the SA method. Usually, the number of absorbers found is small and they can be very 
easily identified in our mock-spectra. However, we checked that the estimated column densities are generally in good 
agreement with the more accurate results obtained with fitting routines. 
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In Figure 6 we show the cumulative number of absorbers per unit redshift range as a function of their column 
density (bottom axis), i.e. the number of systems with column density larger than a given value, and as a function 
of the corresponding equivalent width (top axis), for OVII (left panel) and OVIII (right panel): the continuous lines 
represent the results obtained from hydro-dynamical simulations, while the points are from the SA method. The 
absorption equivalent width, in velocity units, is related to the column density by the simple relation: 

W=j^^XfN (4) 

where N is the column density, / the oscillator strength and A the wavelength of the transition (Saraain 1989, Chen 
et al. 2002). This relation formally holds in the limit of optically thin absorber but is also a good approximation for 
column densities smaller than lO^^cm"'^. 

To quantify the role of the scatter in our simulations, we decided to run another set of SA spectra and switch-off 
any form of scatter, both in the temperature and metallicity. Results without scatter are shown as empty circles, 
while those that include the scatter are represented by filled circles. From Figure 6 we see that, when the scatter is 
considered, the agreement between hydro simulations and SA is good. A first conclusion is that the inputs we used 
for the SA method, taken from hydro-simulations, are enough to have a good description of the number of absorbers 
per unit redshift range for the column density range lO^^'^cm"^ < Novii,oviii < lO^^'^cm"^, corresponding to 
0.1km/s< W< 35 km/s for Ovii- 

The values found are in agreement with the column densities distributions of Chen et al. (2002), which are 
reported in Figure 6 and are represented as a dashed line. In the case in which they take into account the metallicity- 
density relation, parameterized by {logZ/Zq)) = —1.66 + 0.36 log 5 and consider a scatter drawn randomly from a 
lognormal distribution with (log Z/Zq) = —1 and aiogz ~ 0.4, the results are quite similar. We note a slightly 
steeper column density distribution function in our analysis compared to the one in Chen et al. (2002), which results 
in a smaller number of absorbers at intermediate column densities in the range 10^^ — 10^^ cm~^. This is probably 
determined by the fact that our gas-density metallicity relation predicts lower metallicities than the one they used. 

Before making any consideration about the detectability of these systems, we make some comments about the 
effect of the scatter. If we do not take into account the scatter, we cannot reproduce all the systems at column 
densities larger than 10^'' cm"'^ for both the ions. For smaller column densities the effect of the scatter is smaller, 
but it seems that also in this range the scatter produces an increase in the number of systems. 

We checked the different role played by the metallicity and temperature scatter by running another set of SA 
spectra switching off only one of them. It turns out that the scatter that influences most the column density of 
the absorbers is the one in metallicity. This is probably due to the fact that the column density of the absorbers is 
directly proportional to the gas metallicity (Section 2.2). The dependence on the temperature is in a sense less direct, 
as the latter determines the ionization fraction, computed with CLOUDY including both the coUisional ionization 
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Figure 6. The left panel shows the cumulative column density distribution function for OVII (number of absorbers with 
column density larger than a given value per unit redshift). Filled circles are from the semi-analytical model with the scatter 
(i.e. implementing both the scatter in the temperature-density relation and in the density-metallicity relation in the model), 
empty circles are from the semi-analytical model without considering the scatter. The continuous line represents the result 
obtained from hydro-simulations. The dashed line show the results of Chen et al. (2002), for a model which contains scatter in 
the metallicity and a trend of mean metallicity with gas overdensity: (logZ/Z©) = — 1.66 + 0.36 log 5. The right panel shows 
the same quantities for OVIII. Error bars are Poissonian. 



and the photo- ionization contributions (see Section 2.2), which depends on the particular X-ray background chosen 
in a non- trivial way (see for example Figure 4). Another reason could be that the scatter in metallicity is slightly 
larger than the temperature one for the regions with 5 < (5 < 10, a low density WHIM, which is responsible for a 
significant fraction of the absorptions. 

We expect ^ 30 OVII systems with column densities larger than lO^^cm"^ corresponding to equivalent widths 
larger than 2 km/s, per unit redshift range aX z ^ Q. While the number of systems with column densities larger than 
lO^'^'^cm"^, corresponding to equivalent widths larger than ~ 10 km/s, drops to 5. The number of OVIII systems 
is similar to that of OVII ones showing that both ions trace structures at the same density, but the corresponding 
OVIII equivalent width is almost a factor 2 lower. 



4 RESULTS 

To better understand the physical state of the gas which determines the absorptions we plot in Figure 7 and Figure 
8 two AGN pairs, for two different separations of 0.2 h^^ Mpc and 1.2 comoving Mpc (continuous line and 
dashed line represent the two members of the pair) simulated with the semi-analytical technique. Assuming that the 
absorbers are placed at a median redshift of z ~ 0.05, these distances correspond to an angular separation of ~ 3.5 
arcmin and ~ 42 arcmin, respectively, for a ACDM Universe. If the absorbers are placed at a median redshift of 
^; ~ 0.1 the angular separation is ~ 2 arcmin and ~ 30 arcmin, respectively. 

In Figure 7 and 8 it is possible to see that the absorptions arise from gas at a temperature larger than lO'' K 
which has a density larger than lO^'' cm^"^, i.e. overdensity J > 6. One can see that, while for the smaller separation 
(Figure 7) the two AGNs show almost identical features both in density and temperature, which determine very 
similar spectra for OVII and OVIII, for the largest separation of Figure 8 the differences are larger. This is due to 
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Figure 7. From top to bottom: gas density along the line of sight, temperature, spectrum for OVII and spectrum for OVIII. 
Continuous line and dashed line represent the two AGN spectra of a pair with separation rx = 0.2 comoving Mpc at 

z ~ 0.05, simulated with the semi-analytical technique described in the text. This portion of the spectrum has been selected 
around one of the strongest absorptions. 



the fact that now the correlations are weaker, resulting in more different density and temperature profiles for the 
two AGN spectra. Also the spectra now show some coincident absorptions, for example at 5600 km s~^, and some 
anticoincidences, such as the strong one at ~ 10000 km s~^, while at smaller separations there are only coincident 
features (Figure 7). 

To quantify the amount of correlation between the spectra of the AGNs of the simulated pairs we use the 
'hits-and-misses' statistics described by McGill (1990), Crotts et al. (1994), Bechtold et al. (1994), Charlton et al. 
(1997). A coincidence is defined as the case in which an absorption line is present in both the spectra within a given 
velocity difference Au and above some signal-to-noise ratio. An anticoincidence is defined when a line is present in 
one spectrum but not in the other one. If there are two lines within Au, a pretty rare event given the incidence of 
OVII and OVIII absorbers, we count only one coincidence and no anticoincidence, as in Fang et al. (1996). We also 
generate a set of spectra uncorrelated in the transverse direction and compute the same statistics to take into account 
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Figure 8. From top to bottom: gas density along the line of sight, temperature, spectrum for OVII and spectrum for OVIII. 
Continuous line and dashed line represent the two spectra of the AGN pair with a separation of rx = 1-2 comoving Mpc 
at z ~ 0.05, simulated with the semi-analytical technique described in the text. This portion of the spectrum has been selected 
around one of the strongest absorptions. 



the level of random coincidences. Despite the arbitrariness of the definition, which involves the choice of a velocity 
binning for the spectra, it has been shown that it is still useful to define a characteristic size for the observed Lya 
absorbers (Charlton et al. 1997; Fang et al. 1996). The inferred 'characteristic size' of the Lya clouds is in agreement 
with the results of other hydro-dynamical simulations (Theuns et al. 1998) and with semi-analytical models of the 
Lya forest (Viel et al. 2002a, Viel et al. 2002b). 

Results are shown in Figure 9 where we plot the coefficient fco, the fraction of coincidences, which is defined as the 
ratio between the number of 'hits' and the total number of 'hits' and 'misses': the larger the value of this coefficient the 
larger the correlation between the two spectra. These plots have been obtained by simulating absorption spectra of 20 
AGN pairs at a median redshift of z ~ 0.05 for 12 different separations, spanning the range 0.1 — 2 comoving Mpc, 
that is angular separations from 3 to 70 arcmin. OVII and OVIII absorbers are represented by empty and filled cirlces, 
respectively. This coefficient has been computed only for absorbers with column densities > lO^^cm"^, which will be 
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Figure 9. The coefficient representing the fraction of coincidences, as defined in the text, i.e. the number of 'hits' divided by 
the total number of 'liits' and 'misses', plotted as a function of the separation between the two AGNs. Only absorbers with 
column densities larger than 10^"^ cm^^, which can be detected by XEUS and Constellation- X, have been selected to comijutc 
this coefficient. The chosen bin to compute the hits and misses statistics is Aw = 100 km s~^ in the left panel and A.v = 200 
km in the right one. Empty circles represent results for OVII and ffiled ones for OVIII. The dashed line is the level of 
random coincidences, which depends on the velocity bin chosen. The points have been obtained by simulating 20 pairs for each 
separation. Error bars (plotted only for OVII) represent the errors of the mean. 



detectable by XEUS and Constellation- X. Wc have decided to show this result selecting all the absorbers detectable 
by future X-ray missions, i.e. all the absorbers with column densities larger than 10^"* cm^'^. This assumption is not 
crucial: the values assumed by the coefficient fco are pretty robust and change significaiitly only if we include the 
absorbers with lower column densities. In fact, for these absorbers the level of random coincidences is higher and 
determines higher values of fco- The values assumed by fco are mainly determined by the absorbers with column 
densities in the range 

10". 5 _ IQlS 

cm . For larger column densities fco becomes smaller because the number of 

random coincidences found is smaller. 

In the left panel a bin in velocity of 100 km s^^ has been chosen for the hits-and- misses statistics, while in 
the right one the bin is of 200 km s~^. The dashed line is the level of random coincidences calculated from the 
'hits-and-misses' analysis of fiilly uncorrelated AGN spectra. For both absorbers the coefficient peaks within ~ 0.2 
comoving h^^ Mpc and then drops. In both panels the value at which the coefficient becomes comparable to the 
random coincidence level is at ~ 1.2 comoving /i-i Mpc. Of course the number of random coincidences depends on 
the binning chosen, being smaller for smaller A«. However, we stress that the number of random coincidences is very 
small, given the small number of systems expected per unit redshift path. 

Assuming that the noise-level for fco is of the order of 0.05, as we can infer from the right panel of Figure 9, 
a 3a detection of this filamentary structure will require at least 20 pairs with separations ^ 1 comoving Mpc 
corresponding to the range 20 arcmin if the absorbers are placed at zJs 0.1. Another conclusion is that a sample 
of at least 20 AGN pairs with this separation is necessary to discriminate between random coincidences and 'real' 
signal from this filamentary structure. 
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4.1 XEUS and Constellation- X detectability of OVII absorbers 

In this Section we quantify the detectability of OVII absorbers, following the treatment and formalism of Sarazin 
(1989). In particular, we compute detection thresholds of the filamentary network for single AGN spectra and for 
AGN pairs. This moans that in the latter case we will need 20 pairs of sources with separation 20 arcmin, in a 
given portion of the sky, to probe the characteristic size of the absorbers by analysing coincident absorption features 
(see Section 4). We decided to focus only on OVII absorbers. As stressed above OVII and OVIII absorbers probe 
the same structures, but, for a given column density the corresponding OVII absorber equivalent width is a factor 
2 higher than the OVIII one. This means that OVII absorbers can probably be detected in slightly lower density 
regions than OVIII ones, so they can better trace the bulk of the filamentary structure. 

In the following discussion all the fluxes are expressed in the ROSAT All-Sky Survey (RASS) 0.1-2.4 keV energy 
band. 

In order to detect absorption lines, with a specified signal-to-noise level S/N , the detector must collect, within 
the instrument energy resolution AS, a number of continuum photons given by Nphotons^ (■^) ( ^"g" ) j where 
Wabs is the absorber equivalent width. The corresponding flux at the lino energy E is: 

-a ''(ilv)"' {^)\^) (K^)"'(io^)lioob)"'— 

where AQ is the effective collecting area [Q is the quantum efficiency of the detector) and t is the integration time. 

We estimate the source flux needed to detect OVII absorption lines {E = 0.57 kcV, in the observer rest frame) 
using different X-ray telescopes and exposure times of 500 ksec. Although a conservative line detection requires a 
S/N = 5, in all our calculations below we choose S/N = 3. 

We find that the defection of these absorbers is quite difficult. The resolution and collecting areas at ~ 0.5 keV 
of the present generation X-ray satellites are: R = E/ /^E = 500, AQ = 25 cm^ and R = 500 and AQ = 70 cm^ for 
Chandra and XMM-Newton, respectively. If we substitute these values in eq. (5) we find that column densities of lO'^'* 
cm~^, corresponding to equivalent widths Wi^ 10 km s~^, cannot in practice be detected by these satellites: they 
either require exceptionally bright background sources with fluxes ~ 10~®ergcm~^ s~^ (calculated assuming a source 
spectrum F oc E~°''') or unrealistic exposure time of the order of the Chandra life-time. On the other hand, if we 
consider more realistic, although still quite bright and rare. X-ray sources with F ~ lO^^^ergcm"^ s~^ we find that 
the only observable absorption features detectable by Chandra and XMM-Newton are the ones with W > 200 km/s 
(see also Mathur et al. 2002; Fang et al. 2002b). 

Conversely, the low density network of filaments can be detectable by next generation satellites, with much higher 
collecting area, such as Constellation- X * and especially XEUS ^ . For these satellites the effective collecting area and 
resolution are: R = 1500, AQ = 3000 cm^ for Constellation-X and R = 800 and AQ = 40000 cm^ for XEUS. We 
mention here that the value R = 1500 for constcUation-X is that expected from the off-plane option for the Reflection 
Grating Spectrometer (W. Cash 2002, http://constellation.gsfc.nasa.gov) and that for XEUS the ultimate satellite 
configuration XEUS-2 will allow to reach an effective area of AQ = 2 x 10® cm^. 

Assuming that sources with a flux ^ 10~^^erg cm~'^ s~^ arc observed, we find that the minimum detectable 
equivalent width is ~ 6 km s~^ and ~ 3 km s~^, for Constellation-X and XEUS, respectively. With XEUS-2 the 
minimum detectable equivalent width will be ~ 1.5 km s~^. These equivalent widths will probe the range of absorbers 
of column densities of > lO"" cm'^ for Constellation-X, > 10^^-^ cm'^ for XEUS and > 10^^'^ cm'^ for XEUS-2. 

We now estimate the number of AGN with fiuxes sufficient to allow detection of absorbers of given Wabs- We 
use the fit to the log N - log S relation Af{> S) ~ 5.1 x 10~^^(S'/erg cm~^ s~^)~^'®sr~^, originally calculated in the 
Einstein band and here renormalized to the ROSAThand (sec Maccacaro et al. 1982; Hasinger et al. 1993). Formally 
this fit is valid for fiuxes in the range 10"" - 10"" erg cm s and extrapolation to lower fiuxes can be dangerous 
as we know, for example, that the ROSAT Log N - Log S function in the Lockman Hole has a shallower slope than 
the one used here at lower fiuxes (Hasinger et al. 1998). Nevertheless, as we are interested in very bright sources, 
the Maccacaro et al. fit should be a good approximation. Assuming that the typical quasar spectrum is given by 
F oc E~^ '' from eq. (5) and the adopted Log N - Log S relation we find: 

« loHikiv j \—) lioivj (loow^j liooowj liook^j ■ 

In Figure 10 we plot the number of sources per steradian with enough flux to allow the detection of absorbers of 



* http: / /constellation.gsfc.nasa.gov 

t http://astro.estec.esa.nl/SA-general/Projects/XEUS 
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Figure 10. Cumulative source number per steradian with enough flux to allow detection of OVII absorption lines with a 
given equivalent width. The continuous line is for XEUS, the dotted one for XEUS-2 and the dashed one for Constellation-X, 
respectively. We assumed S/N=3 and exposure times of 500 ksec. Vertical lines mark thresholds for which the structures can be 
probed by AGN pairs. The first threshold on the left represents the estimated minimum W detectable for OVII from a sample 
of 20 pairs, with separation ^ 20 arcmin, from an observation of the whole sky. The second threshold represents the same value 
but obtained from a sample of 20 pairs found in a smaller random field of 4 deg^. In the upper axis the corresponding OVII 
column densities are reported. 



a given equivalent width. The solid line refers to XEUS, the dotted one to XEUS-2 and the dashed one on the right 
refers to Constellation-X. In the same Figure we also report the corresponding OVII column densities. 

From Figure 10 one can see that XEUS can detect roughly 10 sources per steradian with equivalent width ~ 2.5 
km/s, corresponding to column densities of 10^'^ * cm~^ for OVII. For these sources, at least a number of 30 features 
can be seen in absorption per unit redshift range if the current models for the WHIM are correct (Figure 8). XEUS-2 
will find a factor of 10 more sources per steradian than XEUS, which will allow to probe the same WHIM structures. 
Constellation-X will not be able to probe the same very low density range. For example, in this case 10 single sources 
per unit solid angle, bright enough to probe absorbers with an equivalent width of ~ 10 km s"'^, are predicted to be 
in the sky. The number of absorption features with this equivalent width is approximately 1 per unit redshift range. 

We estimate now the number of AGN pairs with flux above a given value and whose relative separation is between 
9min and Omax- To this aim we consider a survey area E and we use the X-ray sources angular correlation function 
w{9) = {9/9o)~°-^, with 00 ~ 3 arcmin (Akylas et al. 2000). 

If we indicate with A/e = J^Af{> F) df2 the total number of sources in the survey with fiux above a given value, 
then the number of distinct pairs whose members have a relative separation between 9min and 9max is: 

{l + w{9)) sin 9d9 ^ ^ ^ [Ca. + 1-7 eJ!'** 9l^l,] , (7) 

where the last expression is valid for small separations and for 9min = 0. 

We choose the limits 9min = and 9max ~ 20 arcmin, in computing the number of pairs. For these angular 
separations, as we have found in Section 4, we expect to find significant coincidence features in the spectra of AGN 
pairs. 

The vertical bars in Figure 10 represent the minimum absorber equivalent width which can be probed by using 
AGN pairs, instead of single AGN spectra. Detection of the filamentary structure with future missions can be achieved 
by finding coincident absorptions in AGN pairs, as we showed in Section 4. 

The first threshold on the left in Figure 10 corresponds to the equivalent width which can be probed from a 
sample of 20 pairs, with the right separations, found in the whole sky, while the threshold to the right corresponds 
to a field of 4 deg^ This means that, if XEUS or Constellation-X will observe a random field of 4 deg'^ in the sky, 
they will find 20 AGN pairs with the wanted separation with fluxes high enough to probe the WHIM structures at 
the corresponding equivalent width. If this is the case, the expected number of OVII absorption features seen by 
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XEUS per unit redshift range is ~ 2.5. Thus, if these objects are placed at 2; ~ 0.5 we predict that a sample of 20 
spectra of AGN pairs will possibly show a total of ~ 50 absorption OVII features, while Constellation-X will detect a 
factor of 3 less absorption features. The more powerful XEUS-2 will detect ~ 100 absorptions in the same conditions. 
An observation of the whole sky will possibly allow these satellites to detect WHIM structures at smaller column 
densities. 

To conclude we verify if the RO S AT aX\-sky survey bright source catalogue (Voges et al. 1999) may contain the 
positions of a sufficient number of pairs to probe the WHIM structure. This catalogue contains a sample of 18,811 
sources down to a limiting ROSAT'PSVG count- rate of 0.05 cts in the 0.1 - 2.4 keV band, which corresponds to a 
flux limit 5.2 x 10^^^ org cm^'^ s^^ (assuming a photon spectrum which can be modeled by a power law with slope -1 
and no absorption, 1 cts/s 1 x 10~^^ erg cm~^ s~^). We find that the number of pairs of sources whose separation 
is smaller than 20 arcmin is ~ 100 for fluxes larger than 3.3 x 10~^^ erg cm~^ s~^, and is ~ 300 for fluxes larger 
than 2.2 x lO^^'^ erg cni'^ s^'^. If we consider that approximately 16 % of the X-ray sources have unique counterpart 
of extragalactic origin (Voges et al. 1999), we deduce that the RASS catalogue should contain the positions of ~ 48 
pairs of extrargalactic sources with fluxes ^ 2 x lO"'^^ erg cm ~^ s~^. These bright sources will allow XEUS to probe 
absorbers of ~ 7 km s~^ (see eq. (5)). 

5 DISCUSSION AND CONCLUSIONS 

In this paper wo have presented a semi-analytical technique to simulate X-ray absorption spectra of distant AGNs. The 
Warm Hot Intergalactic Medium (WHIM) can be detected both in absorption and in emission. Recently, a number of 
missions devoted to the study of the WHIM through the observation of emission lines (see e.g. the missions SPIDR 
or IMBOSS ■^) have been proposed. Unlike the lino omission process, which is proportional to the density stjuarc, 
the line absorption has the advantage of being proportional to the density of the intervening absorber. Furthermore, 
absorption measurements are not subjected to background contamination. Absorption measurements, then, happen 
to be more sensitive than emission ones to the lower density regions of the WHIM. Having a larger volume filling 
factor, these low regions constitute the bulk of the filamentary structure. Thus, absorption measurements, besides 
representing an important tool complementary to emission ones, are probably more promising to study the large scale 
structures traced by the WHIM. 

The method proposed here uses as input four results obtained from the z = Q output of a large box-size ACDM 
hydro-dynamical simulations: the probability distribution function of the gas, the gas temperature-density relation, 
the gas metallicity-density relation. Those quantities are sufficient to give a good description of the number of OVII 
and OVIII absorbers expected in the low redshift Universe. Simulations of these two ions are very important: they 
are the strongest trace of the WHIM. If we sum up all the baryons at 2; ~ 2 and compare this value with the baryons 
soon aX z ^ 0, wo discover that a significant fraction of the baryons is missing (Cen & Ostrikor 1999). This fraction is 
thought to be in the WHIM. Thus, the detectability of this network of OVII and OVIII filaments is of fundamental 
importance. 

Wo have extended the simulation of a line-of-sight (LOS) to the simulation of pairs of LOS with the mathematical 
implementation reported in Appendix Al (see also Viel et al. 2002a). This has been done in order to give an estimate 
of the characteristic size of the absorbers by using the information contained in the transverse direction. While 
single LOS arc sensitive to the size of single absorbers, pairs of LOS should bo more reliable a probe of the network 
of filaments traced by the metals as they are more sensitive to the random orientation of each filament. We have 
performed the so called 'hits-and-misses' statistics (McGill 1990; Charlton et al. 1997) on a sample of simulated pairs 
to obtain an estimate of this size. 

The main result of the first part of the paper can be summarized as follows: i) semi-analytical models success- 
fully reproduce the number of absorbers per unit redshift path of hydro-dynamical simulations and we expect ^ 30 
absorbers per unit redshift with an OVII column density > 10^^'^ cm~'^; ii) scatter in the motallicity-gas density 
relation plays a significant role in increasing the number of absorbers with column densities ^ lO^'* cm~^; in) the 
characteristic size of these absorbers at 2 ~ obtained by analysing coincident absorptions in simulated spectra of 
AGN pairs, is ~ 1. 

In the second part we have discussed the detectability of the absorbers by present X-ray missions, such as 
Chandra and XMM-Newton and future ones, such as XEUS and Constellation-X. Up to now only few 'trees' of the 

'X-ray forest' have been detected by Chandra (Fang et al. 2002b, Nicastro et al. 2002a). Most likely, these absorbers 
sample the high column density tail of the distribution shown in Figure 6. A detection of the lower column density 

t http://www.bu.edu/spidr/ 

§ http://www.ias.rm.cnr.it/IMXS/iasimxs.html 
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absorbers, whose volume filling factor is larger, can be achieved only by future X-ray missions whose collecting area 
and sensitivity will be higher. 

We have found that at least a background of 20 pairs of bright extragalactic AGN-likc sources at a median 
redshift of a ~ 0.5 and with relative separations lower than 20 arcmin is needed to probe this filamentary structures 
with 'hits-and-misses' statistics. We have made a theoretical estimate to address the detectability of these sources by 
XEUS and Constellation-X taking into account the bright X-ray sources angular correlation function (Akylas et al. 
2000). The main conclusions are: t) hy observing a random field of 4 deg"^ XEUS and Constellation-X will be able 
to probe WHIM absorbers, with 'hits-and-misses' statistics, which will correspond to OVII column densities 10^*'* 
cm~^ and 10^^'^ cm^^, respectively; ii) in this field XEUS and XEUS-2 will possibly find ~ 50 and ~ 80 absorption 
features, respectively, if we assume that the sources are placed at a median redshift of a ~ 0.5; in) XEUS will detect 
at least 10 sources per steradian bright enough to probe structures with an absorber equivalent width of 2.5 km/s, 
corresponding to an OVII column density of 10^^ ** cm^^; iv) Constellation-X w\\[ detect few thousand sources per 
steradian bright enough to probe structure with an equivalent width of ~ 100 km/s. We have checked in the ROSAT 
all-sky survey bright source catalogue (Voges et al. 1999) to determine if these sources had already been detected. 
We have found that the required number of extragalactic sources with the necessary angular separation is present in 
the catalogue. Re-observation of these sources will allow future X-ray missions to detect baryons in the low redshift 
Universe and to unveil the filamentary structure where they preferentially reside. 
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APPENDIX A: CORRELATION PROCEDURE TO GENERATE MULTIPLE LINES OF SIGHT 
DENSITY AND VELOCITY FIELDS 

We briefly report here the correlation procedure of Viel et al. (2002a). 

Let us start by obtaining auto and cross-spectra for the ID random fields which are defined along single or 
multiple LOS. In the Gaussian case these quantities fully determine the statistical properties of the generated fields. 
If ip{x) is a 3D random field with Fourier transform tpO^) and 3D power spectrum P(|fc|), one can define the LOS 
random field ip\\{x\\,x.±) as the ID Fourier transform 

V>||(fc|||xx) = I l^e*— V(fc||,kx). (Al) 

The cross-spectrum 7r(|A:|| ||rx) for this LOS random field along parallel LOS, separated by a transverse distance r±, 
is defined by 

{V^I|(fe|||xx)V|i(fc|i|xx + rx)) =27rfc(fc|| -hfc||)7r(|fc||||rx) , (A2) 
where 6o is the Dirac delta function and Tv{k\r±) can be related to the 3D power spectrum as follows 

.72 



''^^'^^^ = j • (A3) 

Integrating over angles and shifting the integration variable yields 
1 f°° 

7r(fc|rx) = — dqqMr±y/q^ -k^)P{q) , (A4) 

J k 

where Jo is the Besscl function of order 0. 

In the limit of vanishing distance between the two LOS Jo — » 1 and the above formula reduces to the standard 
relation for the LOS (ID) auto-spectrum in terms of the 3D power spectrum (Lumsden et al. 1989) 
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p{k) = nik\r± = 0) = / dqqP{q) 




(A5) 



Our next problem is how to generate the two random fields and v'G^ in ID Fourier space (—00 < fcy < 00). 
These random fields have non- vanishing cross-correlations but unlike their 3D Fourier space counterparts they cannot 
be related by simple algebraic transformations. 

We gan generally write any M-dimensional Gaussian random vector V with correlation matrix C and components 

dj = (ViVj), as a linear combination of another M-dimensional Gaussian random vector X with diagonal correlation 
matrix, which we can take as the identity I without any loss of generality. The transformation involves the M x M 
matrix A, with components Oij, as follows: V = AX. One gets C = AA^, i.e. dj = ^^,0,^0^^. There is a slight 
complication because V is a random vector field defined in ID Fourier space. Wc can, however, extend the above 
formalism to vector fields, assuming that X is a Gaussian vector field with white-noise power spectrum. 



the conditions X^^.^^ cukajk = Pij- The remaining freedom (due to the symmetry of the original correlation matrix) 
can be used to simplify the calculations. 

It is straightforward to extend our formalism to simulate the IGM properties along parallel LOS. Let V(A:||) and 
W(fc|| ) be two ID Gaussian random vector fields obtained as in Section 2, each with the same set of coefficients Oij 
but starting from two independent whito-noiso vector fields X and Y (i.e. such that (XiYj) = 0). Then both V and 
W have the correct LOS auto-spectra by construction while their mutual cross-spectra vanish: {Vi{k\\)Wj{k'^^)) = 0. 

Let us further define a new vector V'(fc|||rx) with components V/ = X]fe=i (AfcVjt +'yikWk), such that its auto 
and cross-spectra components are given by, 



The vectors V and V' will then represent our physical IGM linear fields on parallel LOS at a distance r±. They 
will be statistically indistinguishable from those obtained by drawing two parallel LOS separated by r± in a 3D 
realization of the linear IGM density and velocity fields. 

The transformation coefficients are determined by the equations "Y^^. i,_-^^{f3ikPji -\- ^ik^jljPki = Pij and 
Sfc=i PikPkj = T^ii- A complete set of coefficients can be found in Viel et al. (2002a). 




J OijXj where 3 of the 4 Oij components are determined by 



(A6) 



mkn)Vj{k\^)) 

(y,(fc||)V;'(fcf||rx)) 



(■K'(fc|||rx)T/'(fc;|jrx)) = 2^ (fc,, + fej, )pi^ 
2tt Soikii + fc||)7rij(|fc||||rx) . 



(|fc||l), 



(A7) 
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